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. ABSTRACT 

We construct a new sample of ~ 1700 solar neighbourhood halo subdwarfs from the Sloan 
Digital Sky Survey, selected using a reduced proper motion diagram. Radial velocities come 
Q from the SDSS spectra and proper motions from the light-motion curve catalogue of Bramich 



^ |et al. I ( |2008| ). Using a photometric parallax relation to estimate distances gives us the full 

C/^ phase-space coordinates. Typical velocity errors are in the range 30 - 50kms"^ This halo 

^ ^ ^ sample is one of the largest constructed to-date and the disc contamination is at a level of 

< 1 per cent. This enables us to calculate the halo velocity dispersion to excellent accuracy. 
We find that the velocity dispersion tensor is aligned in spherical polar coordinates and that 
^ ((Tr, (T0, (T6i) = (143 ± 2, 82 ± 2, 77 ± 2) kms"^ The stellar halo exhibits no net rotation, 

fS| although the distribution of v^ shows tentative evidence for asymmetry. The kinematics are 

consistent with a mildly flattened stellar density falling with distance like r"^-^^. 

Using the full phase- space coordinates, we look for signs of kinematic substructure in the 
stellar halo. We find evidence for four discrete overdensities localised in angular momentum 
and suggest that they may be possible accretion re mnants. The mo st prominent is the solar 
neighbourhood stream previously identified by Helmi et al. ( |1999 ), but the remaining three 



On are new. One of these overdensities is potentially associated with a group of four globular 

O clusters (NGC5466, NGC6934, M2 and M13) and raises the possibility that these could have 

^ been accreted as part of a much larger progenitor. 

•i-H 

Key words: Galaxy: kinematics and dynamics - Galaxy: halo - solar neighbourhood stars: 
^ subdwarfs - Population II 



1 INTRODUCTION 

In recent years, there has been an impressive and growing body of 
evidence that the stellar halo of the Galaxy is composed of the rem- 
nants of accretion and merging events. The Sloan Digital Sky Sur- 
vey (SDSS; [York et al.|2000| ) has made a significant contribution to 
this field, with a series of major discoveries, including the Mono- 
ceros Ring (New berg et al.||2002| [Yanny et al.||2003l ), the Virgo 
Overdensity (|Juric et al. [2008 1 and the Hercules- Aquila Cloud fBe-| 
[lokurov et al.|2007| ). These structures were identified as overdensi- 
ties of resolved stars, exploiting the deep and homogeneous pho- 
tometry in 5 bands (u, g, r, / and z) that SDSS provides in a large 
area around the North Galactic Cap. 

* msniith@ast.cam.ac.uk 



The aim of this paper is to search for substructure in the Galac- 
tic halo with SDSS data, but this time using kinematic methods. 
This technique led to the discovery of the disrupting Sagittarius 
dwarf galaxy and its Stream, which were found serendipitously in 
a radial velocity survey of the outer Galactic Bulge (Ibata et al.| 
|1994| |1995| ). It has also led to the identification of a number of 
probable halo streams found by spectroscopic surveys of the blue 
horizontal branch population (e.g., [Arnold & Gilmore|1992{|Clew-| 
[ley et al.||2005| ). More rarely, searches for kinematic substructure 
have exploited both proper motion and radial velocities. For ex- 
ample, high quality proper motions provided by HIPPARCOS, to- 
gether with ground-based radial velocities, enabled [Helmi et al. [ 
(1999 hereafter H99) to construct three-dimensional velocity dis- 
tributions for an almost complete sample of nearby halo stars and 
to identify a nearby stream as a coherent structure in velocity space. 
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Section Description of cut 



Number of objects 



Density 



1A_ 

n 



Initial sample with high quality astrometry 370,000 

Satisfy RPM cut 27,000 

With high quality spectra 2,704 

Colour within range of photometric parallax 2, 1 83 

Distance < 5.0 kpc 1,717 
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Table 1. The number of subdwarf candidates after the respective cuts. 



Samples of accurate proper motions can be obtained through 
repeated astrometric observations over a significant temporal base- 
line. Most of the SDSS data has limited variability information, 
with one important exception: during the three months when the 
Southern Galactic Cap is available for observation, SDSS repeat- 
edly scanned a ~ 290 square degree area - known as Stripe 82 - to 
detect supemovae (se e e.g.,|Abazajian et al.|200"9|| ). 

|Bramich et al. | ( |2008| ) have presented a Stripe 82 catalogue 
of almost 4 million "light-motion curves", in which objects are 
matched between the ~ 30 epochs, taking into account the eff'ects 
of any proper motion over the eight-year baseline. The catalogue is 
complete down to a magnitude 21.5 in u, g, r and /, and to magni- 
tude 20.5 in z. Each object has its proper motion calculated based 
only on the multi-epoch SDSS J2000 astrometric measurements. It 
reaches almost 2 magnitudes fainter than the SDSS/USNOB cata- 
logue ( [Monet et al.|2003[[Munn et al.|2004| ), making it the deepest 
large-area photometric and astrometric catalogue available. 

Substructure in the catalogue of Bram ich et al. | ( |2008| ) has 
already been identified by looking for overdensities. For example, 
[Watkins et al.| ( |2009j isolated the RR Lyrae variables by a combi- 
nation of colour, metallicity and period cuts. The Hercules-Aquila 
Cloud and the Sagittarius Stream were both clearly identifiable in 
the Stripe 82 RR Lyrae population, together with a completely new 
and distant substructure called the Pisces Overdensity. Thus far, 
however, searches for kinematic substructure in the Stripe 82 data 
have not been carried out, and the object of this paper is to remedy 
this deficiency. 

In Section [2] we show how to isolate a sample of nearby 
halo subdwarfs from the catalogue of [Bramich et al. | ( |2008| ). The 
kinematic properties of the local halo subdwarf population are dis- 
cussed in Section [3] We develop algorithms to search for kinematic 
substructure, recovering the known H99 stream in Section |4] as 
well as new kinematic overdensities in Section [5] 



2 THE SUBDWARF SAMPLE 

2.1 Reduced proper motion (RPM) diagram 



The full [Bramich et al. |p008| ) catalogue contains proper motions, 
jji, for ~ 1 million stars down to magnitudes of r ~ 21.5. We trim 
this sample by enforcing the following three conditions to retain 
only those objects with high quality data. First, we insist that the 
mean object type > 5.7. This requires an object to be classified as 
a star in ^ 90 per cent of epochs. This is less than 100 per cent 
in order to retain objects that have been misclassified in a limited 
number of epochs due to problems with the SDSS star-galaxy sep- 
aration algorithm. In the Bramich et al. (2008 ) catalogue, this oc- 
curs particularly in the final season, when observations were not ex- 
clusively taken in photometric conditions. Secondly, we insist that 
the proper-motion error is less than 4 mas yr~4o remove stars with 
poorly determined proper motions. All of our final sample have 
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Figure 1. Reduced proper motion diagram for the [Bramich et al. [ 12008} 
sample. The solid lines denote the approximate boundary of the subd- 
warf population. The dashed lines denote the stricter boundary employed 
to reduce contamination from disc main-sequence dwarfs (which lie to the 
upper-right of the boundary) or white dwarfs (which lie to the lower-left). 
To aid clarity in this figure, we have imposed an additional proper motion 
cut (// > 30 masyr"^); however, no such cut was applied to our sample. 
Note that the colour-scale saturates at 20 per cent of the peak density. 



proper motions based on at least 23 epochs and time-baselines of 
over 4 years, which indicates that the formal proper motion errors 
should be reliable. Thirdly, we impose r < 19.5 to reduce contami- 
nants (see Appendix A for more details). Here, and elsewhere in the 
paper, all magnitudes are corrected for extinction using the maps of 
Schlegel et al. "('1998V The magnitude cut is often redundant as we 
also later require stars to have SDSS spectra, and there are very few 
spectra for stars fainter than this limit. 

We select our candidate subdwarfs from a reduced proper mo- 
tion (RPM) diagram. The r-band RPM is given by 



Hr-r + 51og- 



mas yr~ 



10, 



(1) 



where ji is the proper-motion and r is the apparent magnitude. This 
uses a star's proper motion as a proxy for distance, allowing us 
to separate cleanly populations with diff'erent absolute magnitudes 
(e.g. main sequence dwarfs, white dwarfs, giants). Although disc 
and halo dwarfs have similar absolute magnitudes, they have very 
diff'erent kinematics. As a consequence, the faster moving halo stars 
appear off'set from the dominant disc stars in the RPM diagram, 
which can be seen clearly in our data (Fig[TJ. 

Note that for the purposes of this figure, we have enforced an 
additional cut that > 30 masyr"^ to emphasise the distinction 
between disc and halo dwarfs. As we relax this criterion, the dia- 
gram becomes populated by slower moving halo dwarfs which lie 
in the boundary region between the two populations on this RPM 
diagram. Slower moving disc dwarfs, on the other hand, lie in the 
upper portion of the diagram well away from our cuts and hence 
do not contaminate our final sample. Therefore, although placing 
such a cut on the magnitude of fi aids the clarity of this figure, it is 
not necessary for our final sample. By not placing a cut on ji, we 
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make the process of quantifying the kinematic bias much easier (as 
described later in Section [24] ). 

We now define our halo subdwarf region as 
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To construct as clean a sample as possible, we reject all objects that 
lie within AH - 0.5 of the boundary, resulting in the final region 
shown by the dashed line in Fig.[T] 

2.2 SDSS spectra 

As well as its photometric survey, the 7th SDSS Data Release 
(DR7; [Abazajian et al.|200 9 ) also has a large number of spectra. 
These are predominantly of galaxies and quasars, but also include 
stellar targets ( [Yanny et al.| [2009 ). These have been analysed by 
a pipeline designed to derive the radial velocities and fundamen- 
tal stellar atmospheric parameters (the SEGUE Spectral Parameter 
Pipeline, or SSPP; ,Lee et al. ,2008a|bl|AlIende Prieto et al . 2008 ). 
Typical internal errors reported by the pipeline are Sy^ ~ 3.8 km s~^ 
for the radial velocity and ^[Fe/H] ~ 0. 1 dex for the metallicity. 
Validation of this pipeline was carried out through comparison with 
high resolution spectra ( Allende Prieto et al.|2008 ), which showed 
that the SSPP external errors are approximately 2.4 km s~^ in ve- 
locity and 0.1 1 dex in metallicity. To account for this, we add these 
errors in quadrature to the internal errors from the SSPP. 

We have cross-matched our candidate subdwarf sample with 
the SSPP DR7 catalogue, and found that ~ 7 per cent have suit- 
able spectra - defined as those with SSPP flag set to 'nnnn', which 
indicates that there are no cautionary signs and the stellar parame- 
ters should be well determined { [Lee et al. |2008a| ). Additionally, we 
reject all spectra with radial velocity error greater than 50kms~^ 
If there are multiple spectra, we take the spectrum with greatest 
signal-to-noise ratio. The median errors in our cross-matched sam- 
ple are 6y, = 4.1 kms"^ and (5[Fe/H] = 0.12 dex. 

The SDSS spectroscopic target selection function is very het- 
erogeneous, covering a wide variety of targets ( [Abazajian et al.| 
[2009[ [Yanny et al.|2009| ). This will introduce selection efl'ects into 
our sample, for example much of the stellar targeting is biased to- 
wards bluer populations which are likely to have lower metallici- 
ties. However, we do not believe that this will introduce any sig- 
nificant kinematic bias. If there was a strong correlation between 
the kinematics and metallicity within the halo subdwarf population, 
then our results would be kinematically biased. However, as we will 
see later in Section [3] we do not find any such strong correlation in 
our sample. 

2.3 Distances, metallicities and surface gravities 

To obtain distances for our subdwarfs, we use a photometric paral- 
lax relation. Such a relation has been established for main-sequence 
stars in the SDSS photometric system by [Ivezic et aL[ ( [200^ . It was 
constructed through observations of eleven globular clusters and 
tested using additional clusters and stars with trigonometric paral- 
lax distances. 

The relation, which is given in equations (A1-A5) of [Ivezic et[ 
[ar] ( |2008l ), is valid for dwarfs with 0.2 < g-i < 4 and incorporates 
a correction to account for systematic trends with metallicity. Ac- 
cording to [Ivezic et al.[ { [2008] ), the intrinsic scatter in this relation 
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Figure 2. Various properties of our subdwarf sample, namely metallicity 
(top left), surface gravity (top right), heliocentric distance (bottom left) and 
spatial distribution in Galactocentric cylindrical coordinates (bottom right). 
Our [Fe/H] distribution may be subject to biases due to the nature of the 
sample construction (see S ection |2. 3). F or comparison, the dashed curve 
shows the distribution from [Ivezic et al. 1 12008 >. The spatial distribution is 
inhomogeneous due to the fact that the proper motion accuracy, which de- 
pends on the number of epochs that a star was observed, varies significantly 
across our field of view (see B ramich et al. |2008| . 



is ~ 0.13 mag, which is a lower limit on the true uncertainty. An 
upper limit of 0.3 has been estimated by [Sesar et al. [ f2008 see 
also Juric et al. 2008), who use candidate binary stars to assess the 
reliability of the photometric parallaxes. Therefore, for this work, 
we take a fiducial uncertainty of 0.2 mag. 

The relation we adopt diff'ers from that of [Ivezic et al.| { [200"8] ) 
in that we have chosen not to incorporate their turn-off" correction 
(given by equation A6 of their paper). Instead, we have constructed 
our own correction using the stellar models of Dotter et al. [ ( [2008] ). 
Our procedure, which is described in Appendix |b[ results in the 
following correction for stars with 03 < g-i < 0.6, 

AMr - aQX + a\xy + a2x' + y + a/^x' y, (3) 

where x = {g-i) - 0.6, y = [Fe/H], and = 2.87, ai = 2.25, a2 = 
-9.79, ^3 = 2.07, ^4 = 0.31. Due to the scatter in this relation we 
increase the uncertainties in the distances of the stars in this colour 
range (see Appendix [B]). 

Using the accurate photometry from the multi-epoch Stripe 
82 data (with mean 6{g - i) < 0.01) and taking [Fe/H] from the 
spectra, we now estimate a distance. Combining the uncertainty in 
the [Ivezic et al.| ( 2008 ) relation with the error on [Fe/H], we find 
that the median distance error in our sample is ~ 1 1 per cent. Our 
distances are mostly in the range 0.5 to 5 kpc, although there are a 
small number of stars further out. However, unless otherwise stated, 
henceforth we only use stars with heliocentric distance D < 5 kpc. 
Beyond this, the errors on the velocities become too large because 
of the uncertainty in the proper motion. This gives us a final sample 
of 1717 halo subdwarfs, as summarized in Table [T] 

There are a number of potential sources of contamination in 
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Figure 4. Multi-component Gaussian fit to the velocity distribution. The 
statistically preferred fit is this two-component model, which is better able 
to reproduce the asymmetry than the single Gaussian shown in Fig. [3] 



the sample, including white dwarfs, thin and thick disc stars and 
background giants. These are discussed in turn in Appendix A, 
which concludes that the actual level of contamination is very small 
(< 1 per cent). As a check that our sample is clean, we examine var- 
ious non-kinematic properties in Fig. |2] We can immediately see 
that the surface gravities are consistent with that of a dwarf popu- 
lation. Although there are a handful of stars with low gravity, the 
vast majority have log(g) > 3. These outliers are most likely caused 
by misestimation of log(g) from the spectra. The median error on 
log(g) for our sample is 0.25 dex, but the median error for those 
stars with log(g)<3 is 0.5 dex. 

The distribution of [Fe/H] is also consistent with a halo sam- 
ple. In Fig. |2] we show the observed distribution along with the 
Gaussian found by jlvezic et al.| (2008 ) using photometric metallic- 
ities. Although the mean of our distribution ([Fe/H] = -1.55) is 
similar to theirs ([Fe/H] = -1.46), it is clear that our distribution 
is significantly broader (with a dispersion 0.43 as opposed to 0.3). 
This broadening and offset is most likely due to the underlying se- 
lection function for the SEGUE spectra, which will clearly have a 
significant effect on the [Fe/H] distribution. 



2.4 Kinematic bias 

One problem with our sample is that it suffers from kinematical 
bias. This is caused by the cut on H^, which selects stars via their 
tangential velocity, rather than their proper motion. Using a cut on 
//r, instead of the proper motion, makes the task of quantifying the 
bias significantly easier, since assumptions about the underlying 
distance distribution (i.e., luminosity function) are not needed. 

We calculate our detection efficiency as follows. For each sub- 
dwarf in our final sample, we take the sky coordinates and create 
a mock sample of 50,000 fake stars. Then, for each mock star, we 
select Mr and g-/ at random from our observed distribution. Note 
that for each realization both the magnitude and color are assigned 
from one star. We then fix the kinematics by choosing velocities 



from Gaussians with the halo velocity dispersions fromjKepley et| 
[aT7] < |2007| ) and no net motion (folio wing p^llende Prieto et al.|2006l L 
We determine the tangential component of the velocity (vtan) with 
respect to the line-of- sight between the Sun and the mock star. Q 
Given this information we can calculate the reduced proper motion 
using. 



(4) 



The efficiency is then given by the fraction of mock stars which 
pass our reduced proper motion cut. 



3 GLOBAL KINEMATIC PROPERTIES 
3.1 Data 

Here, we calculate properties of the first and second moments of 
the halo velocity distribution. We begin by defining our coordinate 
systems. We take the solar radius as 8 kpc, the velocity of the lo- 
cal standard of rest as 220 km s~^ and the solar peculiar velocity as 
given by pehnen & Binne3r| ( |1998] ). The Galactocentric Cartesian 
reference frame is denoted by (x, y, z), where the axes are oriented 
along the line connecting the Sun and the Galactic centre, in the di- 
rection of disc rotation and toward the North Galactic Pole, respec- 
tively. This is a right handed frame so that the Sun is at x = -8.0 
kpc. Velocity components resolved with respect to this coordinate 
system are (Vx, Vy, v^). 

We also use cylindrical and spherical polar coordinates de- 
fined with respect to a right-handed Galactocentric frame. Cylin- 
drical polars are denoted by (R, 0, z) where R is radially outward, 
is positive in the direction of counter-rotation of the disc and z 
is positive towards the North Galactic Pole. The corresponding ve- 
locity components are (vr, v^, v^). Similarly, spherical polars are 
denoted (r, 6, 0) where r is radially outward and 6 is increasing to- 
wards the South Galactic Pole. The velocity components are (Vr, v^, 
v^). So, for stars in the solar neighbourhood, v^ ^ -v^. Note that 
disc stars rotate with v^ ^ -220kms~\ which is a consequence of 
adopting a right-handed system. 

The values for the mean and dispersion of the velocity compo- 
nents in spherical and cylindrical polars are given in Table [2] They 
are calculated using a maximum likelihood technique which cor- 
rects for the spread that results from observational errors, namely. 



1 



exp 



(5) 



where fit and cr^ are the mean and dispersion of velocity compo- 
nent /, whilst Vi^k and Svi^k are the velocity and its uncertainty for 
the ^h star, 6k is the detection efficiency for this star, and N is the 
total number of stars in the sample. The errors on the individual 
velocities are calculated using Monte Carlo methods, incorporating 
the errors on the proper motion, distance and radial velocity. The 
median error in (Vr,V0,Ve) is (38, 47, 35) kms~^ The magnitude 
of the errors differs between the three components because of the 
location of our field, which causes the proper motion to contribute 



To calculate Vtan requires us to assign a distance to each mock star, which 
we do at random from the observed distribution. This means that the effi- 
ciency does have a dependence on the distance distribution; however this 
dependence is negligible since equation J4| is a function of the tangential 
velocity rather than the proper motion. 
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System 


Component 


<v> (kms-i) 


cTy (kms ^) 


Spherical 


Vr 

ye 


39+2.8 
2.3+1I 
^- -1.4 


142.7+2-0 
82.4 ± 1.4 
77.3 ± 1.1 


Cylindrical 


Vr 

Vz 


8.9 ±2.6 
-1.2+1.6 


138.2+2-0 
89.3+|-3 



Table 2. Efficiency weighted mean velocities and their dispersions for our halo sample in both spherical and cylindrical polar coordinates. The values are 
obtained using a standard maximum likelihood technique which incorporates the spread that results from observational errors (see equation [5}. The cross- 
terms or covariances are reported in a separate paper (Smith, Evans & An |2009) . 



Reference 



System Component <v> (kms cry (kms ^) 



WooUey 



1978 



Spherical 



ye 



148 
122 
82 



Norris 



1986 



Cartesian 



131 ±6 
106 ±6 
85 + 4 



Gould & Popowski |{T998 i 
(Kinematically selected sample) 



Cartesian 



-5± 15 
9± 13 
1 ± 8 



171 ± 11 
99 ±8 
90 + 7 



Gould & Popowski 



1998 



(Non-kinematically selected sample) 



Cartesian 



-3 ±9 
28 ±9 
2±5 



160 + 7 
109 + 9 
94 + 5 



Chiba & Beers 



1 2000 



Cartesian 



-16 ± 16 
26+ 12 
-5± 11 



141 ± 11 
106 ±9 
94 ±8 



Gould 



j2003a 



Cartesian 



-1 ±3 
9 (fixed) 
2±3 



171 ± 10 
99 + 8 
90 ±7 



Kepley et al. ] ( |2007: 



Cylindrical 



Vr 



-4± 11 
-23 ±8 
-1 ±6 



157 ±8 
110 + 6 
84 ±4 



Table 3. Some determinations of the halo velocity parameters from the literature. Note that v^ is defined counter to the direction of disc rotation and for nearby 
samples v^ ^ -Vy. With the exception of Kepley et al. ( 2007), the estimates for the bulk motion have been standardised using the value of the solar motion 
from |Dehnen & Binney ||T998) . [Kepley et al. H2007| use the value from |Mihalas & Binney ||l98T) and as a consequence their estimate of (v^) is subject to 
an offset of ~ 7kms~\ 



Component 


Corr[v/, Vj] 




[Vr, V^] 


0.052 ± 0.028 


2.2 ± 1.2 


[Vr,V0] 


-0.019 + 0.042 


-1.6 ±3.6 


[v^,ve] 


-0.073 + 0.050 


-36.7 ±24.1 



Table 4. Efficiency weighted correlation coefficients and tilt angles for halo stars with | z | > 1 kpc (see |Smith, Evans & An |2009| for the definitions of these 
quantities). Note that these values are not identical to those quoted in Smith, Evans & An [(2009) because the ones given here incorporate our new turn-off" 
correction for the photometric parallax relation (see equation [3}. However, the diff^erence is small and none of the conclusions of [Smith, Evans & Air|(2009t 
are aff^ected. 
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Figure 3. The velocity distributions of the halo subdwarf sample. The velocity components are resolved in spherical polar coordinates with prograde rotation 
corresponding to v^ < 0. The solid curves show Gaussian fits to the data, with the corresponding parameters given in Tablel2| 



more to the component. Furthermore, because the uncertainties 
are larger in the tangential direction the errors are correlated, with 
the orientation of the error ellipse varying along the field. 

The maximum likelihood velocity distributions are displayed 
in Fig. |3] Although the Gaussian models provide a reasonable 
match to the observations for Vr and Nq, it appears that v^ exhibits 
some asymmetry. Accordingly, we repeated the maximum likeli- 
hood fitting allowing for multiple Gaussian components. We deter- 
mine whether a multiple Gaussian model is statistically preferred 
by comparing the ratio of likelihoods, where -2 ln(L//L/+i) > 3 
indicates that a model with (/ + 1) components is preferred. The 
value of 3 comes from the fact that there are three additional pa- 
rameters in the model with (/ + 1) components. 

Although for Vr there is only a marginal improvement with 
multiple Gaussians, for v^ a two-component fit is clearly preferred, 
reflecting the fact that the distribution is narrower than a single 
Gaussian yet has noticeable wings for \wo\ > 200kms~^ As a con- 
sequence, this two-component model has one main narrow compo- 
nent (with cr = 64kms~^) and a shallow broad component (with 
cr = 135kms~^), both of which have mean close to kms~^ The 
presence of these wings can be understood when one considers the 
substructures present in our sample (which will be discussed later). 

For V0, a two-component fit is also preferred, with approxi- 
mately equal contributions from the two components with disper- 
sions of 48.1 and 98.9 kms~^ and means 22.8 and -17.1 kms~\ 
respectively (see Fig.|4]). This provides a better match to the asym- 
metry in the v^ distribution. We believe that this asymmetry is not 
caused by problems with our efficiency correction nor with contam- 
ination from the disc (see Appendix A). This leads us to conclude 
that it is a real effect. It appears that the asymmetry is more pro- 
nounced for metal-rich stars ([Fe/H] > -1.5 dex), although we 
find no evidence for dramatic trends with metallicity, such as those 
postulated by [Morrison et al.| ( |2009| ). Neither do we see any clear 
gradient in (v^) as a function of height from the Galactic plane. Al- 
though we identify various kinematic substructures in our sample 
(see Sections|4]and[5]), none of these can explain this asymmetry. 

In Table [3] we present earlier determinations of the mean ve- 
locity and dispersion for halo stars in the literature. Our dispersions 
are significantly smaller than previous estimates (such as Gould 
Popowski |1998| [Chiba & Beers |2000l [Gould ^20Q3a, Kepley et, 
al. |2007| ), although the ratios are in approximate agreement. Com- 
parison between different investigations is clouded by the fact that 
some of the earlier samples are subject to significant levels of disc 
contamination. Also, it is not always clear whether corrections have 



been made for the spread induced by measurement errors. The dis- 
persions estimated by |Gould"| { [2003a ) are upper limits as they do 
not incorporate the spread due to uncertainties in their photometric 
parallax relation. 

There has also been controversy in the literature with regards 
to the mean rotational velocity of the halo, with estimates of pro- 
grade (Chiba & Beers 2000, Kepley et al. 2007 ), no rotation 
( Gould & P opowski 199 8 , Allende Priet o et al. 2006) , or retro- 
grade (Maj ewski [1992 1. Our value of (v^) = 2.3^}-3kms"^ is ap- 
proximately consistent with a non-rotating halo. Since our kine- 
matic selection is biased against stars with v^ ^ -200kms~\ we 
might expect that our estimate of (v^) is dependent on the efficiency 
correction. However, tests show that this is not the case. Our cal- 
culation in Section [2.4[ requires an assumption for the rotational 
velocity of the halo, but if we adopt values anywhere between -20 
and 20 kms~\ then our measured value of (v^) varies by less than 
1 km s~^ However, it is still true that our estimate of (v^) (as with 
any of the estimates from Table [3]) is degenerate with the assumed 
value for the local standard of rest (220 km s"^. 

Another result of interest is the value of (Vr), which is 3cr away 
from zero. It is not clear why we obtain a positive value of (Vr) 
(i.e., radially outward). There are a variety of potential issues that 
could cause this effect: kinematic substructure; systematic errors 
in the velocities (either the radial velocities or proper motions) or 
distances; the presence of binary stars, for which we would system- 
atically underestimate their distances; bias in the determination of 
the solar motion with respect to the local standard of rest. We do 
not investigate this issue further in the present paper. 

We also calculate the efficiency weighted correlation coeffi- 
cients and tilt angles, as defined in [Smith, Evans & An [ ( [2009) , for 
the 1,568 stars with |z| > 1 kpc. Although the data used here is 
the same as that of [Smith, Evans & An [ ( [2009] ), our sample diff'ers 
in that we are now incorporating our new turn-off correction to the 
parallax relation (see equation |3] and Appendix [B]). The values are 
given in Table |4] These are very similar to those quoted in [Smith] 
[Evans & An (2009 ) and their conclusions are not affected by this 
new parallax correction term. 



3.2 Modelling 

Taking our cue from the recent arguments of [Fellhauer et al.[ ( [2006| ), 
let us assume that the gravitational potential of the dark halo is 
spherically symmetric. If the halo has a flat rotation curve of am- 
plitude vo (= 220 kms~^), then the potential is 
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i/r = -v^logr. 



(6) 



We now seek a phase space distribution function for the stellar halo 
that can reproduce the observed kinematics of the SDSS subdwarfs. 
From Jeans Theorem (see, e.g., Binney & Tremai ne |1987| l, the dis- 
tribution function must depend on the integrals of motion, namely 
the binding energy E, the components of angular momentum, J^, 
Jy, Jz, together with the total angular momentum, /. 



A number of authors { |White|1985l|de Zeeuw et al.|1996| ) have 
shown that the stellar density laws of the form 



p = por ^ sm 

can be reproduced by distribution functions 

f^,,(E, J\ Jl) = 7/^,,72'"/f exp[(r+2f7+2m)£]. 



(7) 



(8) 



where m+n > -1 and 2n> -\, and rfm^n is a normalization constant 
(given in equation 3.5 of |de Zeeuw et al.|1996) . Note that we have 
written m and n so that our notation is consistent with earlier work, 
but m and n are not necessarily integers. 

The corresponding velocity second moments are 



P(vl) 



2m + 2n + y 
m + n + \ 



rJ 



n ■ 



1 



(9) 



pivj) = (2n + l)p{vl). 



So, for example, if m = ;i = 0, the model is isotropic and all three 
velocity dispersions are just vq/ ^fy. This is a well-known result 
for isothermal populations. More generally, for fixed y and n, the 
velocity dispersions and the density have the same angular depen- 
dence for all m, but the anisotropy ratios (v^>/(v^> and (v^>/(v^> do 
depend on m. 

A flexible way to model the stellar halo is to build further dis- 
tribution functions by linear superposition. Here, our aim is to con- 
struct a simple distribution function that reproduces the kinematics 
of the SDSS subdwarfs, and so we choose 



f(E, J ,J ) = O^m,0fm,0 + (^m,lfm,U 



(10) 



where the ani,n are constants specifying the fraction contributed by 
each component. The corresponding density law is flattened and 
has the form 



p{r,e): 

We require that at = 8.70 kpc and z ■ 
our SDSS subdwarf population) 



(11) 

-2.41 kpc (the centroid of 



(vi) 



.2\l/2 . 



143kms 



.2\l/2 . 



: 77kms 



(vi) 



.2\l/2 



82kms-\(12) 



The advantage of a two component model ([To| is that there is a 
unique solution. There are three velocity dispersion constraints, and 
there are three unknowns, namely the radial fall-off" y, the index m 
and the ratio pi/po (or equivalently CKm,i/Q^m,o)- Solving the non- 
linear simultaneous equations numerically gives the solution y = 
3.75, m = -0.72 andpi/po = 0.063. 

In other words, the SDSS subdwarf kinematics are consistent 
with a stellar halo in which the density falls off" like r~^^^, somewhat 
steeper than the r~^^ advocated in the classical work on the metal- 
poor populations of the halo (e.g. , ,Freeman| 1987 The axis ratio of 
the stellar distribution can be computed from 



1 + 



P]_ 
Po 



-l/y 



(13) 
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Figure 5. The distribution of angular momentum for our subdwarfs (left) 
and for an artificial smooth halo (right). These plots cover the distance range 
- 2.5 kpc (top) and 2.5 - 5 kpc (bottom). The contours denote our detec- 
tion efficiency in steps of 10 per cent, where the solid line corresponds to the 
90 per cent contour. The box corresponds to the location of the kinematic 
stream identified by |Helmi et al. | ( p999) . In this coordinate system the Sun 
would lie at approximately (-1800, 0) kpc kms~^. 
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Figure 6. The distribution of angular momentum residuals for subdwarfs in 
the distance range - 2.5 kpc (top) and 2.5 - 5 kpc (bottom). The residuals 
are obtained after subtracting a smooth model from the observed distribu- 
tion. The contours and H99 boxes are as in Fig. [5] In this coordinate system 
the Sun would lie at approximately (-1800, 0) kpc kms~^. 



which gives q = 0.983, in other words, very round. This distribu- 
tion function is not unique, as there are undoubtedly more com- 
plicated multi-component models. Nonetheless, it is the simplest 
distribution function that is consistent with the kinematical data. It 
is interesting that the triaxial kinematics are consistent with a near- 
spherical stellar density law. 
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Figure 7. The distribution of the angular momentum components of our 
subdwarfs in the distance range - 2.5 kpc (top) and 2.5 - 5 kpc (bottom). 
Stars are colour-coded according to [Fe/H] and the colour-range saturates at 
-2.5 and -0.5. Contours denote the residuals shown in Fig.|6] corresponding 
to 1, 3, 5 and 10 sigma significance. There are a number of stars which lie 
outside the range of this plot; these are given in Table |5] In this coordinate 
system the Sun would lie at approximately (-1800, 0) kpc kms~^ 



4 KINEMATIC SUBSTRUCTURE 
4.1 Quantification of Substructure 

The components of the angular momentum perpendicular and par- 
allel to the symmetry axis of the Galactic plane are 



[(yVz-zVy)H (zvx-xvz)^l' 



Jz = xvy-yvx. (14) 



For a solar neighbourhood sample, these two angular momentum 
components are essentially ~ -roVy and ~ r© |Vz|. Our median 
errors on (Jz, J±) are (410, 255) km^ s~^ If we restrict our sample 
to the 645 subdwarfs within 2.5 kpc, then the errors are reduced by 
a factor of ~ 35 per cent. Use of these coordinates is a common 
approach (e.g. H99, |Chiba & Bee rs 2000), as they are adiabatic 
invariants in a spherical potential. Note that, unlike H99, we use a 
right-handed coordinate system, so our values of take an oppo- 
site sign to theirs. This means that in our coordinate system the Sun 
would lie at approximately (-1800, 0) kpc km s~^ . We weight each 
star according to the inverse of its detection efficiency to account 
for kinematic bias and smooth the distribution using a Gaussian 
kernel with FWHM of (300, 150) kpc km s~\ which is the same or- 
der of magnitude as the errors. The resulting distribution is plotted 
in Fig.|5] 

To test for substructure, we generate an artificial sample of 
halo stars using the observed distances and directions of stars in our 
sample, but generate velocities according to the trivariate Gaussian 
with means and dispersions given in Table [2] We smooth this dis- 
tribution and subtract it from the smoothed observed distribution 
(see Fig. [5]). The resulting residual plot is shown in Fig. [6] where 
the colour-scale shows the significance of overdensities. The sig- 
nificance is quantified by generating a series of efficiency-corrected 
realisations from our distributions, all of which have the same num- 



ber of stars as the observed sample. The scatter in these realisations 
provides an estimate of cr(Jz,J^), from which we can deduce the 
significance of our overdensities. We show the overdensities for 
two distance ranges (the 645 stars with D < 2.5 kpc and the 1072 
with 2.5 < D < 5 kpc) and include contours showing the detection 
efficiency. 

Some of the apparent overdensities in the outer regions are 
artefacts. The smooth model predicts very few stars in the outer re- 
gions and hence a-(Jz,J^) is close to zero, which exaggerates the 
significance of individual outliers. For example, the large clump 
around (Jz, J^) = (2000, 2500) kpc km s~^ in the upper panel of Fig. 
|6]is actually composed of only one star. This can be seen on com- 
parison with Fig.|7] which shows the unbinned (Jz, J±) distribution 
for each of the stars. Note that the figures confirm that our sample 
is free from significant disc contamination. If there were significant 
thin or thick disc stars, they would be visible in the region around 
-2500 < 7z < -1500 kpc kms-^ and < 600 kpc kms-^ 

4.2 The H99 Kinematic Stream 

The most noticeable feature in Fig. [T] for stars within 2.5 kpc is 
the asymmetry in Jz for stars with 1500 < J^ < 2500 kpc kms~^ 
There is a significant number of stars located around (Jz,J±) ~ 
(-1000, 2000) kpc km s~^ . This corresponds to the solar neighbour- 
hood kinematic stream first identified by Helm i et al. |(T99 9). They 
analysed a sample of nearby (i.e. < 1 kpc) halo stars with full six- 
dimensional phase space information and found that ~ 10 per cent 
come from a single coherent structure. Subsequently, [Kepley et al.| 
||( |2007 ?. hereafter K07) probed a larger volume {D < 2.5 kpc) and 
found a smaller fraction (~ 5 per cent) of stars belonging to this 
stream. Restricting their sample to 1 kpc, they find a larger fraction 
(~ 9 per cent) in agreement with H99, although this is not sur- 
prising since there is significant overlap between the H99 and K07 
samples. Chiba & Beers (2000) also studied an extended version 
of the H99 halo sample and found a smaller fraction than H99, al- 
though the exact percentage is unclear. K07 defined the location of 
the H99 region as -1500 < Jz < -500 and 1400 < J^ < 2500, 
which is indicated in Figs. [5] -[T] 

There has also been related work on streams in the solar neigh- 
bourhood subdwarfs by Gould (2003b). This analysis, which was 
based on a large sample of objects with proper motion measure- 
ments but without radial velocities or distance determinations, used 
a Bayesian likelihood analysis to investigate various properties of 
the nearby subdwarf population. Instead of looking for direct evi- 
dence of cold kinematic streams, he quantified the amount of gran- 
ularity in the velocity space and concluded that no more than 5 per 
cent of local halo stars (within ~ 300 pc) can come from a single 
cold stream. Although this number appears to be in conflict with 
the original results from H99, there is actually no discrepancy as 
the H99 stream is formed from two separate clumps in velocity 
space. 

The H99 stream is clearly visible in our data, as can be seen 
in the top panel of Fig.|6] although it appears more diffuse than the 
detection in K07 (c.f., their fig. 11). Note that this difference cannot 
be explained by the size of our observational errors, since the me- 
dian error on (Jz, J^) for our H99 stars is (240, 209) km^ s~^ . Closer 
inspection of the upper panel of Fig.[7]shows that the H99 box con- 
tains two stars near the centre along with a more diff'use component 
of 10 stars. These two stars lie almost exactly in the centre of the 
box at (Jz, J±) ~ (-1035, 1954) kpc kms~^ The mean metallicity 
of the two stars is -2.2 dex, which is more metal poor than our full 
subdwarf population (-1.5 dex). However, the mean metallicity of 
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Figure 8. Properties of our H99 stars (see Section [4!2) . Upper panels and 
lower-left panel show the velocity distributions for all stars within 2.5 kpc 
(small dots) and stars in the H99 region (large circles). The bottom-right 
panel shows [Fe/H] as a function of offset from the centre of the H99 box 
in the 2D (J^, J^) plane, which lies at (-1000, 1950) kpc km s"^ . 



all 12 stars in the H99 region is -1.8 dex with a dispersion of 0.4 
dex, which is comparable to the overall distribution. 

The stars in and around the H99 region are illustrated in Fig. [8] 
which shows their velocities and offset in the (J^, J^) plane from the 
centre of the H99 box (the mean error in Joffset is 196 kpc kms"^. 
Given that a smooth halo distribution should give ~ 4 stars in this 
region, we suspect that not all of these 12 stars are associated with 
the H99 feature. The distance distribution of these stars is approxi- 
mately uniform between ~ 1 and ~ 2.5 kpc, with no obvious trend 
towards smaller distances (unlike K07, who found 67 per cent of 
their H99 stars were within 1 kpc). However, we note that our sam- 
ple is biased against D < 1 kpc, as can be seen from Fig. [2] which 
implies that we could be missing stars from this feature. If we con- 
sider the more distant stars with D > 2.5 kpc as shown in the lower 
panel of Fig. [7] we find further potential members of the H99 re- 
gion. However, these stars appear to come from a more homoge- 
neous background of stars with high angular momentum. 

We wish to quantify the prevalence of the H99 stream. We first 
need to correct for the detection efficiency in order to account for 
the kinematic bias in our sample. However, the efficiency in Section 
|2.4| is not directly applicable for stars in the H99 stream, as they do 
not have the same velocity distribution as the underlying smooth 
halo. To account for this we re-calculate the efficiency for stars with 
H99 kinematics. Using this efficiency we predict that up to ~ 7 
per cent of all halo stars within 2.5 kpc lie in the H99 box, which 
is significantly larger than the fraction expected from an entirely 
smooth halo (~ 1 per cent). 



4.3 Retrograde and Prograde Outliers 

Various authors have noted the presence of prograde or retrograde 
outliers when analysing distributions of halo velocities. For exam- 
ple, j^nnetaT] ([20^ identified halo stars on extreme retrograde 



orbits (V0 > 200kms~^) and found them to be significantly more 
metal poor than typical halo stars. K07 similarly saw an excess 
of retrograde stars in their halo sample and noted a deficiency in 
[Fe/H] (they found 6 stars with ([Fe/H]> = -2.1 dex). K07 argued 
that these stars probably belong to a tidally disrupted stream. In ad- 
dition, K07 found 3 metal-poor stars on extreme prograde orbits, 
although they postulated that these were more consistent with the 
underlying halo distribution. 

In our sample, we find many stars on extreme orbits, as can be 
seen from Fig. |7] There are also a number of stars whose angular 
momenta is sufficiently large for them to lie outside the region of 
this plot; those stars are listed in Table [5] 

However, none of the low outliers in Fig |7] look particu- 
larly clumped and neither the prograde nor retrograde outliers are 
significantly in excess of expectations of a smooth distribution. 
If we analyse stars with < 2000, then we find 5 stars with 
< -2000 kpc kms"^ and 31 stars with > 2000 kpc kms'^ 
Our smooth model predicts 2.4 and 30.3 stars in these regions, re- 
spectively. Furthermore, neither of these two samples are partic- 
ularly metal-poor; the mean metallicity is -1.6 dex and -1.4 dex 
for the prograde and retrograde outliers, respectively. Note that the 
strong asymmetry in the number of outliers is due to the fact that 
our detection efficiency is significantly lower for < 0. 



5 NEW KINEMATIC OVERDENSITIES 

There are a number of other possible clumps in our sample. We dis- 
cuss three of them, which we label Sloan Kinematic Overdensities 
(SKO), in the following section. The associated stars are given in 



Table[6] 



5.1 SKOa 

One interesting feature in the nearby sample shown in the upper- 
panel of Fig. |7] is the collection of 3 retrograde stars centred on 
(Jz,J±) ^ (410, 2420) kpc km s"^ The three stars are all more 
metal-poor than the halo average, with ([Fe/H]> of -2.1 dex, and 
are not clumped spatially (see Table [6]). The stars are enclosed by 
a circular region of the (Jz,J±) plane centred on (414.8, 2420.3) 
with radius 350 kpc km s~^ . Our smooth model predicts that there 
should be only 0. 1 per cent of the sample in this region, which is 
significantly smaller than the observed fraction of 0.5 per cent. We 
shall refer to this potential overdensity as SKOa. 

When constructing our subdwarf sample, we enforced a strict 
cut on the quality flags that are raised by the SSPP (see Section [2!2| . 
The SSPP includes two categories of flags: critical and cautionary. 
For the latter category, it is often possible to determine the stellar 
parameters, although they should be treated with caution. The vast 
majority of these stars are flagged because the SSPP pipeline raised 
issues regarding the Ha line strength ( |Lee et al. [2008 a| ). In partic- 
ular, this flag marks stars where 'there exists a strong mismatch 
between the strength of the predicted Ha line index, based on the 
measured US line index' . If we include stars with cautionary flags 
into our subdwarf sample, then we gain an additional ~ 700 stars, 
of which three appear to be associated with the SKOa feature. All 
of these stars are have similar metallicities to the SKOa members 
mentioned above (see Table [6|. 

We plot the velocities of the six SKOa stars in Fig. [9] along 
with the distribution of [Fe/H] as a function of heliocentric dis- 
tance. Also included in this figure are the stars in a looser selection 
region with radius 500 kpc km s~^ and with distances up to 5 kpc. 
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15 
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-1 
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17 
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14 
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11 
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12 
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-1 
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12 
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Table 5. Retrograde and Prograde Outliers 



SKO SDSS ID Jz Jx Vr V0 \e Distance [Fe/H] 
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a 
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-63.2 ± 42.8 


99.9 ± 102.0 


-245.5 ± 86.1 


6 


31 


± 





85 


-1 


17 


± 





11 


bt 


SDSS J030948.29-000111.3 


1371.6 ± 1650.7 


2433.0 ± 1455.8 


15.6 ± 49.7 


102.7 ± 121.0 


-162.2 ± 108.3 


7 


87 


± 





95 


-1 


61 


± 





12 


bt 


SDSS J03 1120.73-000329.4 


1655.1 ± 1097.4 


2735.3 ± 1050.5 


245.6 ±37.1 


141.2 ± 89.6 


-216.2 ± 76.9 


5 


45 


± 





70 


-1 


10 


± 





16 


c 


SDSS J015840.78-000115.6 


1305.3 ± 690.3 


2002.9 ±461.6 


-177.6 ± 38.9 


130.4 ±67.5 


-183.8 ±47.4 


4 


12 


± 





40 


-2 


61 


± 





17 


c 


SDSS J015840.90+002900.8 


1367.9 ± 864.8 


2889.3 ± 729.4 


87.6 ± 49.5 


133.1 ± 82.3 


-258.5 ± 62.3 


4 


60 


± 





52 


-2 


61 


± 





18 


c* 


SDSS J015806.00+003419.5 


1825.3 ± 1378.6 


2688.8 ± 1247.4 


228.8 ± 62.3 


162.2 ± 119.1 


-201.9 ±97.3 


6 


55 


± 





82 


-2 


51 


± 





12 



Table 6. Stars which are candidate members of the three new kinematic overdensities (SKOa-c). Note that the errors on the velocities and angular momenta 
can be highly correlated; they are included here to give an indication of the relative uncertainties. Notes: * The Sloan Stellar Parameter Pipeline flagged the 
parameters for these stars as 'cautionary' due to issues regarding the Ha line strength; ^ The distances for these stars are greater than 5 kpc, which means that 
the derived velocities are subject to large uncertainties; ^ The Sloan Stellar Parameter Pipeline flagged the parameters for this star as 'cautionary' due to the 
presence of a strong G-band feature in the spectrum, and the distance is greater than 5 kpc and so the velocities are uncertain. 



It seems that a number of the more distant stars may also be asso- 
ciated with this overdensity. 

Although the above discussion of SKOa appears tentative, 
there is striking evidence supporting an accretion scenario for these 
stars. In Fig.[T2] we overplot the locations of globular clusters using 
data from Dinescu, Gir ard & van Altena ,(19 99 ) on our distribution 
of kinematic overdensities in the (Jz, J^) plane. We immediately see 
the clump of four clusters standing out around SKOa. These clus- 
ters are NGC5466, NGC6934, NGC7089/M2 and NGC6205/M13, 
which have metallicity -2.22, -1.54, -1.58 and -1.65 dex, respec- 
tively. Of these four clusters, NGC5466 is known to be disrupting 
( [Odenkirchen & Grebel |2004l|Belokurov et al.|2006l ) and its orbit 
takes it to within a few kpc of the solar neighbourhood (Fellhauer, 
private communication; see also Fellhauer et al. ,2007 ). This, along 
with the fact that the SKOa stars have similarly low metallicity, 
indicates that they could well be tidal debris from NGC5466. How- 
ever, the fact that there are four clusters located in the same region 
of the (Jz, J^) plane raises the intriguing possibility that all could be 



associated to a single major accretion event. Previous authors have 
questioned whether subsets of these four clusters may be associ- 
ated (e.g. [D inescu, Girard & van Altena '1 999[[Palma, Majewski &] 
[Johnston ^2 002 , Mackey & Gilmore 2004|), but it appears that no 
one to date has attempted to explicitly associate all four clusters. 

The uncertainties on the velocities of these clusters varies; 
NGC6205 and NGC7089 have reasonably well constrained angu- 
lar momenta (with error on Jz and J^ of between 200 and 400 
kpc kms"^), while for NGC5466 and NGC6934 this is less well 
constrained (with error of between 400 and 900 kpc kms"^. This 
is due to the uncertainties in the proper motion combined with the 
large distances to these two clusters (~7 to ~15 kpc). However, 
although the error for NGC5466 seems large, its proper motion 
is confirmed by Fellhauer et al. [ p007[ ) and so its location in the 
(Jz, J±) plane is robust. 



Note that a similar retrograde feature has been identified by 



Dinescu 


d2002 


see also Brook et al. 2003) from the sample of 


Chiba & Beers 


([2000b. This feature, which they postulated could 
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Figure 10. Properties of the SKOb overdensity (see Section js^ and Table[6}. The coloured symbols denote stars with > 2000 kpc kms~^ and are colour- 
coded according to metallicity, with SKOb members (large symbols) differentiated from other high-J^ stars (small symbols). The triangles denote stars with 
D > 5 kpc, which have larger uncertainties in their kinematics. We also show stars with < 2000 kpc kms~^ for comparison (small points). The vertical 
dashed lines denote the right ascension range 46° < a < 48° . Note the lack of high-J^ stars outside of this range. 
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Figure 11. Properties of SKOc sample stars in the range 28° < a < 30° (see Section [5^ and Table[6j. The upper-middle and upper-left panels also show stars 
outside this a range for comparison. The coloured symbols denote stars with > 2000 kpc km s^ and are colour-coded according to metallicity, with the 
two metal-poor SKOc members (large symbols) differentiated from other high-J^ stars (small symbols). 
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Figure 9. Properties of the SKOa overdensity (see Section [sT] and Table 
|6j. The upper two panels and lower left panel show the velocity distribu- 
tions for all stars within 2.5 kpc of the Sun (small points), members of the 
SKOa overdensity (large circles), additional members with cautionary SSPP 
flags (triangles), and potential members with distance greater than 2.5 kpc 
(small circles). The members with D < 2.5 kpc are enclosed within a 350 
kpc kms~^ circle centred on (413.6, 2419.9) kpc kms~^ while the poten- 
tial members with 2.5 < D < 5 kpc are enclosed within a larger radius (500 
kpc kms~^). The bottom right panel shows [Fe/H] as a function of heho- 
centric distance for stars with ofl'set in the 2D (Jz, J±) plane less than 500 
kpc kms~^ from the centre of the SKOa clump. 
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Figure 12. Angular momentum distribution of the 38 globular clusters from 
[Dinescu, Girard & van Altena (1999) with 6D phase-space, where the 
colour corresponds to cluster metallicity. Contours denote the kinematic 
overdensities for halo stars within 2.5 kpc (see Section|4j, corresponding to 
1, 3, 5 and 10 sigma significance. Note the four globular clusters potentially 
associated to the overdensity at (Jz, J±) ~ (500, 2500) kpc km s~^ ; these are 
NGC5466, NGC6934, NGC7089/M2, NGC6205/M13. In this coordinate 
system the Sun would lie at approximately (-1800, 0) kpc kms~^. 



velocities of |Chiba & Beers | ( |2000| ) are significantly smaller than 
ours and no more than one of their stars lies in our SKOa region. 



5.2 SKOb, SKOc 

It is clear from Fig. [T] that there are a significant number of more 
distant stars (2.5 < D < 5 kpc) with large values of J^. To es- 
tablish whether these stars are clumped, we must investigate their 
spatial distribution along the stripe. As our sampling is very inho- 
mogeneous, we are actually interested in identifying locations on 
the stripe where the fraction of high objects is anomalous, rather 
than the absolute number. 

Since our field of view consists of a narrow 2.5° wide stripe 
at constant declination, streams will typically cut through the stripe 
and be localised in right ascension a. Therefore, we analyse our 
data by looking at 2° sections in a along the stripe, concentrating 
only on bins which contain more than 30 subdwarfs. We calculate 
the fraction of stars in a bin with > 2000 kpc kms~^ (consid- 
ering all stars with D < 5 kpc), and also the expected fraction as 
predicted from our smooth model. When we carry out the calcula- 
tion of the binomial probability, we take into account the fact that 
there are 21 bins with more than 30 stars, i.e. we determine the 
probability that any of these 21 bins has the observed fraction of 
high stars. The smooth model predicts a factor of ~ 1.9 fewer 
high stars in total compared to the observed fraction (60/1717). 
This indicates the presence of non-Gaussian tails in the velocity 
distributions and it is these outliers that could harbour potential ac- 
cretion remnants. 

There are two locations on the stripe where the binomial prob- 
ability of obtaining the observed fraction of high stars is less than 
~0.05 and for which there are more than three such stars. We call 
these two clumps SKOb and SKOc. Candidate members of these 
overdensities are given in Table |6] 

The clump SKOb is located in the range 46° < a < 48°, where 
we find 8 of 57 stars have > 2000 kpc kms~^ Our smooth 
model predicts that we should see only 2.1 stars and the corre- 
sponding binomial probability is 0.03. Various properties of the 
SKOb stars are shown in Fig. [To| from which we see that there 
are a significant number of stars in this a range at distances be- 
tween 3 and 5 kpc. It is intriguing to note that most of these objects 
have retrograde velocities (Jz > 0). There is a clear clump of three 
stars located in a small region of the (v^, v^) plan^ around (89, 
-225) kms~^ Furthermore, these three stars have similar metal- 
licities (([Fe/H]> = -1.5 dex) and distances (see Table [6|. Given 
that these stars are close to our limiting distance of 5 kpc, it is 
worth looking at those stars which lie outside this distance cut. 
So in Fig. [T0| we also include stars in this a range with distances 
up to 8 kpc. Although their velocities are less certain, there are 
sufficient stars to warrant their inclusion. Clearly, these more dis- 
tant stars reinforce the significance of the potential overdensity. 
The clump in the (v^, v^) plane is even more pronounced, with six 
stars now lying in this region. All six stars are located in a loose 
clump in the (Jz, J±) plane, with Jz e (670, 1660) kpc kms~^ and 
J^ G (2430, 3250) kpc km s"^ Note that these additional three stars 
also have similar metallicities to those with D < 5 kpc. If we loosen 



be debris torn from the system that once contained oj Centauri, 
was observed as an excess of stars with small retrograde velocities. 
However, as can be seen by comparing our Fig.[9]to fig. 15 of |Chiba| 
|& Beers | ( [2000| ), the stars in our SKOa feature are not compatible. 
Although our stars indeed have small retrograde velocities, the J^ 



^ Note that if these stars were debris from an accreted sateUite, one might 
naively expect them to also have the same values of Vr. However, as can 
be seen from the modelling of satellite disruption for the H99 stream, this 
is not necessarily the case; although the model presented in Kepley e t al. | 
j2007| is localised in (Vz, V0)-space, it has an extended vr distribution. 
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the cut on the SSPP flag as above (see Section [5T| ), then we do not 
gain any additional potential members for this overdensity. 

The second clump, SKOc, is located at 28° < a < 30° and 
has 7 of 69 stars with > 2000 kpc kms~^ The smooth model 
predicts only 1.9 stars (with corresponding binomial probability of 
0.07). The properties of these stars are illustrated in Fig.[TT] Unlike 
SKOb, there is no obvious clumping in velocity space. However, 
of note is the pair of two metal-poor stars which lie at very sim- 
ilar distances (D ^ 4.1,4.6 kpc, [Fe/H] ^ -2.6,-2.6 dex) and 
are relatively close to each other on the sky (separated by half 
a degree). These two stars are located in close proximity in the 
(V0, v^) plane, but their angular momenta are not particularly close 
(J^ ^ 2000 + 460,2890 + 730 kpc kms'^. This means that the 
chance of these two stars belonging to a kinematic overdensity is 
perhaps unlikely. Despite this, their case is strengthened to some 
extent when we include stars with less secure parameters; there is 
one additional metal-poor star located in the (v^, v^) plane at (162, 
-202) kms"^ with D ^ 6.6 kpc and [Fe/H] ^ -2.5 dex. As well 
as being more distant than our standard cut of 5 kpc, there is also 
a cautionary flag raised in the SSPP. The flag records that the star 
may be exhibiting a strong CH G-band (around 4300 A) relative 
to what is expected for a 'normal' star ( |Lee et al. |2008a| ) and im- 
plies that the spectral parameters should be treated with caution. It 
lies very close on the sky to the pair of metal-poor stars previously 
mentioned, but given the various uncertainties it is far from clear 
whether we can claim that it is associated to the pair of metal-poor 
stars in this region. Therefore, in conclusion, the validity of this 
feature is rather tentative and it could simply illustrate that we have 
reached the limits of this method. 



6 CONCLUSIONS 

This paper has presented a catalogue of ~ 1700 halo subdwarfs 
within a heliocentric distance of ~ 5 kpc. Our analysis is restricted 
to a particular patch of sky. Stripe 82 of the Sloan Digital Sky 
Survey (SDSS), and exploits the light-motion curve catalogue con- 
structed by lBramich et al. |(2008| ). The exceptional precision of the 
SDSS photometry allows us to determine distances with an uncer- 
tainty of ~ 10 per cent from a photometric parallax relation. As 
a consequence, we obtain high-quality tangential velocities (with 
typical errors of 30 - 50kms~^) using the proper motions derived 
from the light-motion curves. Radial velocities are obtained from 
SDSS spectra. So, this sample of halo stars has full six-dimensional 
phase-space information, allowing us to probe in detail the kine- 
matical structure of the stellar halo. As our data are restricted to 
Stripe 82, we may be susceptible to spatial variations present in the 
halo population. However, even if the halo population is not well- 
mixed, our ~250 square degree field-of-view is sufficiently large to 
overcome any small-scale eff'ects. 

The size of our sample allows us to determine the orientation 
and semi-axes of the velocity ellipsoid of the stellar halo to ex- 
cellent accuracy. The velocity ellipsoid is almost precisely aligned 
in the spherical polar coordinate system, indicating that the total 
potential must be nearly spherical ( [Smith, Evans & An ||2009| ). 
The velocity dispersions are (cr^, cr^, cr^) = (143 + 2, 82 + 2, 77 
+ 2) kms~^ The values of the dispersions are significantly smaller 
than previous estimates, although the ratios of the dispersions are 
in good agreement. A simple distribution function that matches the 
kinematic data is constructed - and it suggests that the density of 
the stellar halo is only mildly flattened and falls ofl" with distance 
like p ~ r~^-^. 



The stellar halo exhibits no net rotation. The velocity disper- 
sions are reasonably well fit by a Gaussians, with the exception of 
the V0 component which is asymmetric. We believe that this is a 
real eff'ect. It appears to be more pronounced for metal rich stars in 
our sample. None of the kinematic substructures are able to account 
for this asymmetry. 

On the issue of substructure, we confirm the presence of the 
stream identified by |Helmi et al. |(T99 9 ). Using the angular momen- 
tum components (Jz, J±), we find a number of other potential sub- 
structures, which we label Sloan Kinematic Overdensities (SKOs). 
Two of these new features are particularly striking, containing a 
number of stars that are localised in both kinematics and metal- 
licity. The metal-poor overdensity SKOa (([Fe/H]) ^ -2.1 dex) 
appears to be coincident in the (J^, J^) plane with an association of 
four globular clusters (NGC5466, NGC6934, NGC7089/M2 and 
NGC6205/M13), suggesting that they may have been part of the 
same accretion event. If so, then this implies that the progenitor 
must have been a large satellite, similar in size to Fornax. We are 
currently investigating the orbits of these stars and clusters and plan 
to test this hypothesis using simulations. 
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APPENDIX A: POTENTIAL CONTAMINANTS 

There are three main sources of contamination: white dwarfs, disc 
dwarfs and background giants. We discuss each in turn and con- 
clude that our contamination level is less than ~ 1 per cent. 

Al White dwarfs 

In a RPM diagram, the white dwarfs are fainter and bluer than the 
halo subdwarfs. Their presence in Fig.[T]is only just discernible as 
the cloud around ((g-i),Hr) ~ (-0.4, 17). Although the density of 
white dwarfs in this figure is clearly small, it is possible that white 
dwarfs with small proper motions enter our subdwarf region. As 
stated in Section[2^ we require that the SSPP flag is set to 'nnnn'. 
This flag includes a check for candidate white dwarfs f Lee et al7] 
2008a ) and so our sample should be free from such contamination. 

As a check, we cross-matched our subdwarf sample with the 
catalogue of spectroscopically confirmed white dwarfs from the 
fourth data release (DR4) of the SDSS ( [Eisenstein eFaTII^OOgl ). 
This catalogue was constructed via both automated criteria and vi- 
sual inspection and so should be reasonably complete. We find that 
none of their white dwarfs overlap with our subdwarf sample. Al- 
though the |Eisenstein et al. | ( |2006| ) catalogue contains only DR4 
data, we know that ~ 60 per cent of our subdwarf spectra were in- 
cluded in DR4 and so, assuming the fraction of white dwarfs does 
not vary with SDSS data releases, we can expect < 1 white dwarf 
to contaminate our subdwarf sample. 

A2 Disc dwarfs 

Although our RPM cut in Section[2T| rejects stars with kinematics 
similar to the Sun (hence retaining mainly fast moving halo stars), 
it is possible that some disc stars could enter our sample. Contam- 
ination by thin-disc stars is negligible due to two factors. First, the 
small scale-height (~ 200 pc) of the thin disc means that our sam- 
ple, which goes out to a few kpc below the plane, should have min- 
imal contamination. Second, the thin disc is cold so that it is highly 
unlikely that a thin disc star would have large enough proper mo- 
tion to pass our RPM cut. However, it is possible that there may 
be some contamination from the hotter thick-disc component. Pre- 
vious studies based on similar RPM cuts have found levels of disc 
contamination to be of order two per cent ( [Gould |2007| ). 

To quantify the contamination from disc dwarfs in our sam- 
ple we analyse the kinematics of the stars. We calculate the prob- 
ability that a star belong to a particular population (i.e. thin disc, 
thick disc or halo) by comparing its velocity to a toy model con- 
structed using Gaussian distributions for each component of the 
velocity, taking into account the diff'erent scale-heights and local 
density normalisations. For the scale-heights and disc kinematics 
we follow |Smith et al. | ( [2007| ), with the exception of the mean rota- 
tion velocity of the thick disc which we model using the parameters 



Kinematics ofSDSS subdwarfs 15 



determined from |Girard | ( |2006| ), i.e. = -195 + 30|z| kms'^ We 
model the halo component using the kinematics of |Kepley et al. | 
pOOV) and assume that the halo density does not vary significantly 
within the volume probed. As expected, we find that none of our 
stars are compatible with the thin disc. For the thick disc we anal- 
yse the stars with [Fe/H] > -1.5 dex, since this is the three-sigma 
lower limit on the thick-disc metallicity according to | Soubiran"etl 
p. (2003 ). Although there are no stars in this [Fe/H] range with 
thick-disc membership probabilities of greater than 0.36, there are 
six stars with probabilities of between 0.2 and 0.36. If we add up the 
thick-disc probabilities of all stars in this [Fe/H] range we can es- 
timate the total level of contamination. Note that this will probably 
be an overestimate because the kinematics of the thick disc overlap 
those of the halo, which implies that even a pure halo sample will 
be estimated to have a non-zero level of contamination. Using this 
technique we estimate that the thick-disc contamination should be 
no more than 10.2 stars, i.e. less than one per cent. This should have 
little bearing on our results. As a test we repeated the analysis of 
Section [3] using only stars at distances of more than 2.5 kpc from 
the Galactic plane and found that global kinematic properties were 
unchanged. 
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A3 Background Giants 

Another possible source of contamination is from faint background 
giants. However, given the RPM cut, we see that a star must have 
a sufficiently high 11 to shift it into the subdwarf regime. For exam- 
ple, a stationary distant star with r - \9 mag and a spurious ji - A 
mas yr~^ would have Hr - \2 mag. If the star was blue enough then 
it could conceivably pass the Hr cut. However, in our final subd- 
warf sample we find that the star with the least significant proper 
motion is inconsistent with = at the 2.5cr level, which means 
contamination from distant (and hence stationary) stars is negligi- 
ble. 



APPENDIX B: TURN-OFF CORRECTION FOR 
PHOTOMETRIC PARALLAX RELATION 

We estimate distances using a photometric parallax relation based 
on [Ivezic et"aL] ( |2008| ), who construct a relation using data from 
a number of globular clusters. They first determine a (g - /) - My 
colour-magnitude sequence for stars on the main- sequence by iden- 
tifying the colour of the main- sequence turn-off for these clusters 
(i.e. the point at which the slope of the colour-magnitude relation 
becomes vertical) and discarding all data within 0.05 mag of the 
turn-off. Once they have this relation, they then devise a correction 
to account for the presence of the main- sequence turn-off, which 
they base on the sequence for the cluster Ml 3. 

However, the morphology of the colour-magnitude relation 
around the turn-off region will clearly be dependent on both metal- 
licity and age, and so we would like to construct a correction which 
incorporates these effects. In order to do this, we use the stellar 
models of [Potter et al. [ ( [2008] ), taking isochrones with ages from 1 
to 15 Gyr (in steps of 0.5 Gyr) and [Fe/H] from -0.5 to -2.5 (in steps 
of 0.5 dex). Following the compilation of [Venn et al. | ( |2007] )m we 
choose [of/Fe] =0.3, which is representative for our halo [Fe/H] 
range. 

We take each set of isochrones and, in the same manner as 
[Ivezic et aT] pOOSj ), we shift each sequence so that it has = 
for {g-i) = 0.6. Then, for a given [Fe/H], we calculate mean Mr as 
a function of (g-i) in the range 0.3 < (g-i) < 0.6, considering only 



Figure Bl. Our correction to the photometric parallax relation due to the 
main-sequence turn-off. The points denote mean magnitudes from the stel- 
lar models and the corresponding solid curves denote the polynomial fit 
given in equation [B7] with [Fe/H] increasing from -2.5 (bottom line) to -0.5 
(top Hne). The dashed Hne corresponds to the turn-off correction of |Ivezic| 
|ernH2008) . 

model data up to the main- sequence turn-off. Since we know the 
approximate age distribution of the halo population, we calculate 
a weighted mean using a Gaussian prior with mean and sigma of 
10 and 2 Gyr, respectively.!^ We then calculate the offset between 
the mean model magnitude and the uncorrected relation of |Ivezic| 
|et al. 72008'^, i.e. that given in equations (A1-A5). This is shown in 
Fig. |Bl| where the effect of the turn-off correction is to increase the 
brightness of stars compared to the uncorrected relation. 

This data is then fit with a third-order polynomial, resulting 
in the following equation which is applicable for stars in the range 
0.3 < (g-i) < 0.6 and -2.5 < [Fe/H] < -0.5, 

AM^^ - X + ai xy + a2 x' + a^x^ y + xy^, (Bl) 

where x = (g-i) - 0.6, y = [Fe/H], and ao = 2.87, ai = 2.25, ^2 = 
-9.79, ^3 = 2.07, ua - 0.31. Note that we do not include terms 
depending solely on y since we require the relation to go through 
AMJ^ = for (g-i) = 0.6. We also discard the x^ term since this 
has no effect on the fit. Although the fit is far from perfect, given the 
overall uncertainties in the method we believe this should provide 
a reasonable approximation. 

In order to avoid any divergent behaviour in this relation we do 
not extrapolate equation beyond [Fe/H] = (-2.5, -0.5). For 
stars with metallicities outside this range, we use the relation at the 
limit of our allowed range (i.e. -2.5 dex for the metal-poor stars and 
-0.5 dex for the metal-rich stars). This has very little bearing on our 
results: although we have around 30 stars with 0.3 < (g-i) < 0.6 and 
[Fe/H] < -2.5, the relation converges at the metal-poor end and so 
the [Fe/H] = -2.5 relation should provide a good approximation; 

^ One could also incorporate a prior based on the initial mass function or, 
equivalently, the luminosity function. However, we investigated this and 
found it had a negligible effect on our results. 
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from Fig. |Bl| it can be seen that the relation does not converge at 
the metal-rich end, but this is of little consequence since there are 
fewer than 10 stars with 0.3 < (g-i) < 0.6 and [Fe/H] > -0.5. 

Clearly the uncertainties on the parallax relation will be larger 
in this colour range due to the scatter in this turn-off correction. 
To estimate the uncertainty, we calculate the standard deviation in 
Mr when we calculate the mean. The scatter varies as a function 
of colour and metallicity, but if we take the relation with the largest 
scatter (corresponding to [Fe/H] = -0.5 dex) we find the following 
relation, 

S [amJ^) = 0.39 - 0.65 (g-i). (B2) 

When estimating distances for our stars we add this uncertainty in 
quadrature to the other sources of error. 



